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An eikonal expansion is used to provide systematic corrections to the eikonal approximation 
through order where k is the wave number. Electron wave functions are obtained for the Dirac 

equation with a Coulomb potential. They are used to investigate distorted-wave matrix elements 
for quasi-elastic electron scattering from a nucleus. A form of effective-momentum approximation 
is obtained using trajectory-dependent eikonal phases and focusing factors. Fixing the Coulomb 
. . . distortion effects at the center of the nucleus, the often-used ema approximation is recovered. Com- 

OO \ parisons of these approximations are made with full calculations using the electron eikonal wave 

, functions. The ema results are found to agree well with the full calculations. 
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I. INTRODUCTION 



Professor Manoj Banerjee was a collaborator and friend of both of the authors. It is an honor to contribute a paper 
that is dedicated to his memory. He exhibited an enthusiasm and passion for physics that has inspired many. 

In this paper, we summarize some methods and results concerning the issue of Coulomb corrections in quasi-elastic 
scattering of electrons by nuclei. Electrons have been used extensively as an experimental probe of the internal 
structure of nuclei in the past few decades. Electron scattering is considered to be a precise tool in view of the 
strength of the e.m. interaction, so that as a result the scattering process can be treated in the one photon exchange 
^ ' approximation. In particular, quasi-elastic scattering of electrons by nuclei has been used to investigate properties 
Ch ! like the validity of the Coulomb sum rule in nuclei. Exp eriments have been performed at the MIT Bates Laboratory 

' ^ ■ 0, i, i, i, i, i, SHU, at the Saclay Laboratory [11 [O, E El, H El and at SLAC [M E3, El ■ Although the 

I \ electron may in general be considered a weak probe, complications arise due to Coulomb distortion effects in the 
^ 1 electron wave function owing to the nuclear charge distribution. As a result, in order to extract nuclear and nucleon 
\^ ' structure information from these experiments, the Coulomb distortion contributions have to be accounted for in the 
0^ , theoretical analysis of the data. In addition there exists the complication of the presence of final state interaction. 
' Studies have shown that in the quasi-free region at high momentum transfer these effects are expected to be small. 
, Neglecting the final state interaction, the Coulomb distortion can be handled in the so-called distorted wave Born 
approximation. Here the Coulomb potential is treated exactly by solving numerically the Dirac equation in the 
presence of the Coulonib potential. Exact solutions for the Dirac-Coulomb wave functions may be obtained as a sum 

■ over partial waves. [l9l. [20l. [2l| . As the electron energy increases, the partial- wave expansions converge more slowly in 

■ spite of the fact that the Coulomb corrections become smaller. These calculations are numerically complex and have 
^ . the disadvantage of not allowing for a simple theoretical interpretation in terms of nuclear structure functions. So 

' one important theoretical issue has been to investigate whether there is a simple yet reliable way to characterize the 
reaction process in terms of response functions, similar to what can be done in the absence of final state interactions and 

■ Coulomb distortions. Confining attention to high enough energy, a natural and reliable framework for the description 
is given by the use of the eikonal wave function. 

The eikonal approximation clearly gets better at increasing energy and it allows a simpler analysis of the effects of 
Coulomb corrections. Some particulary transparent results have been obtained using the eikonal approximation to 
derive an effective-momentum approximation (ema) p^ . [2I that produces results very similar to plane- wave results. 
It is important to include focusing factors such as those found in the WKB approximation [23| and revisited more 
recently in quasi-elastic scattering ^26. , 27.] . However, the attempts to combine the eikonal analysis with focusing 
factors suffer from the lack of a systematic basis. Significant disagreements in the determination of nuclear response 
functions from experimental data (28l . [2^ have arisen at least in part owing to the use of different theoretical methods 
to remove the Coulomb corrections. Therefore, it is of interest to study a systematic expansion of the eikonal approach, 
where the various effects arise in a natural way. 

In order to address the issue of Coulomb corrections, we have developed corrections to the eikonal approximation 
based upon a systematic expansion in the high-energy limit (30| . This eikonal expansion is shown to be rapidly con- 
vergent already at typical energies of few hundred MeV for targets used in quasi-elastic scattering. Corrections to the 
eikonal approximation have a long history. Work by Saxon and Schiff [31] showed how to correct the approximation 
to leading order in A systematic expansion for the scattering t-matrix was developed by Sugar and Blankenbe- 
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cler [s^ . Systematic corrections to the Glauber approximation [33*1 were developed in [s^, [s^ and extended to the 
Dirac scattering amplitude in |,3Q] . However, a systematic expansion for wave functions has not been developed prior 
to Ref [3^. 

In Section |TT1 we present the eikonal expansion for the Dirac wave function and show that the focusing effect can 
be obtained at order \/k of the expansion. In particular, we focus on u(r), which is a Pauli spinor containing the two 
upper components of the Dirac wave function. The lower components are simply ±1 times the upper components 
because of helicity conservation. Convergence of the eikonal expansion is shown to be fast in the few hundred MeV 
electron energy region. Because there is a spin-orbit interaction, spin-dependent terms arise in the eikonal expansion. 
They also are determined and their effects are found to be negligibly small. In Section IIIII we summarize the basic 
formulae for quasi-elastic electron scattering. In Section IIVI we deal with effective-momentum approximations and 
discuss the original ema approach. A natural modification (EMAr) to the ema approximation is proposed, where 
trajectory-dependent eikonal phases and focusing factors are included. 

Specializing to the longitudinal response, section |V] discusses quasi-elastic scattering by use of a simple model of 
the nuclear response. Comparisons of the full calculations of the response functions are made with the two effective- 
momentum approximations. In the fits to the plane- wave impulse approximation with an effective momentum small 
deviations from unity of the normalization are found. However, overall reasonable agreement is found with the full 
calculations, lending support to the use of these effective-momentum approximations as a basis for the theoretical 
analysis of the quasi-elastic data. Some concluding remarks are made in Section [VII 



II. THE DIRAC ELECTRON WAVE FUNCTION 



In this paper we consider electron scattering at intermediate energies, where the eikonal approach is expected to be 
reasonably accurate. The electron is assumed to be a Dirac particle. For the Dirac equation, the eikonal expansion 
is carried out in two stages. First we consider the Pauli spinor ■u(r) that contains the two upper components of the 
Dirac wave function, i.e.. 



V(r) 



u{r) 



It follows from the Dirac equation that the Pauli spinor £{r) that contains the two lower components may be determined 
in a second stage, where the two lower components are determined in terms of u(r). 

Eliminating the lower component spinor from the Dirac equation we find for the upper-component spinor the 
equation 



^E,-V,-a-p^—y^a-pju{r), (I) 

where Ei = E — m, E2 — E + m and E is the energy of the incoming particle. For electron scattering it is generally 
the case that E » m and thus Ei k, E2 ^ E. Because the electron mass is much smaller than the energy, helicity 
is conserved and the lower components are given simply by ^A(r) = 2\u\{v), where A = ±1/2 is the helicity. 

For outgoing-wave boundary conditions, the Pauli spinor u{y) is written in terms of a complex eikonal phase 
x'^' = x'^n^") + i'^'^^H^) E^nd a complex spin-dependent phase 7*^+^ = 7(+'(r) + i(5^+-'(r) as follows 



(r)=(l--^) f,^I^^^^X^•'f,^<y.T•\ (2) 



where k — \J E^ — -m? is the momentum of the incoming wave. The wave propagates in the z-direction and an impact 
vector b is defined as the part of r that is perpendicular to the z-direction, i.e., b = z x (r x z). Three orthogonal unit 
vectors are : z,b = b/|b| and e = bx z. The spin matrix in the eikonal phase is — a ■ e. The factor (1 — Vc/E2y^'^ 
is introduced in order to sum up terms that otherwise arise in higher orders. 

The eikonal expansion has been developed in Ref. [3^. The result is that the eikonal phases are expanded in a 
systematic fashion in powers of 1/k as 

xW = x(+^+xi+^+xr^ + --- 



,M = 7W+7(+) 



= Si+^ + .--, (3) 
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where the subscript of each term denotes the power of \/k that is involved. ExpUcit expressions can be found in 
Ref. |3fl]. The leading terms are given by 



(+). ^ 1 r , idVcir) 



2k db 

S^+\r) = 0, (4) 

where higher order terms than 1/fc have been dropped. 

The upper-component spinor of the Dirac wave function for helicity A and outgoing-wave boundary conditions is 
given by 

4+) = //^W(r)e^'=^e^x*^'e'-=^'^'6, (5) 
where is a helicity eigenstate. The Dirac focusing factor J^^^^ is defined as 

/f^"^(r)=(l-|-yV.'^'. (6) 

One may work at various orders of the eikonal expansion by truncating the expansions of Eq. ([3]). Similarly, the upper- 
component spinor for helicity A and incoming-wave boundary conditions can be obtained by replacing in Eqs. (j4l6p 
the superscripts by (— ) and the integration ranges 

/-co by ST- 

Convergence of the eikonal expansion has been studied for scattering of a 500 MeV electron. A rough estimate can 
be made of the higher order corrections of the expansion. Given that the electron mass is m = .511 MeV, it follows 
that k K, E and w « 1, both within a part per million. The Coulomb potential is approximately Vc(0) — 25 MeV at the 
center of the nucleus. The eikonal expansion introduces corrections that involve the nondimensional ratio Vc/ E w .05, 

so we expect ^^^^ « .0025 It should be noted, that the eikonal expansion is not convergent but is asymptotic, meaning 

that the error should be bounded by the first neglected term. 

Figure [T] shows the eikonal phases for a charge Z = 100 and electron energy E — 200 MeV. The Coulomb potential 
is chosen to be 

where — Vo is the value of the potential at r = and i? is a range parameter. This Coulombic potential corresponds 
to a charge density 

_ 3^0 

^^"■^ 47re (r2 -f i?2)5/2 ' 

The above parameters were chosen in order to make the corrections visible. The corrections are much smaller for a 
500 MeV electron and a smaller nuclear charge. 



III. QUASI-ELASTIC ELECTRON SCATTERING 

Let us consider the quasi-elastic nucleon knock-out process {e,e',N) from a nucleus. The cross section of quasi- 
elastic electron scattering can be expressed in terms of the transition matrix element Ai : 

''dn.^WEjpE,, (9) 



d^fdEf J (2n) 



where p is the momentum of the knocked-out nucleon and Ep = ^JM^^^^^-pp- is its energy. The bar denotes an average 
over initial hclicitics and a sum over final helicities. We have 



M = / [ I J,^(q, p) (10) 



I I , \ I I I I I \ , \ I I 

1 2 3 4 5 6 7 
b(F) 



FIG. 1: Eikonal phases at z=0 versus impact parameter: solid line shows xo, dashed line shows xi, dotted line shows X2, 
rectangles show cui and ovals show uj2. A constant has been added to xo such that it vanishes at 6 = 3.5R. Phases are shown 
for Z=100, E=200 MeV and R=2 fermi. 



with the electron current matrix element 

= Jd\ui-^*(r)j,e-^^-'^ui+\r), (11) 

for emission of a photon of energy lo = Ei — Ef and momentum q = ki — kf . In Eq. (fTTj) u^x^~^ are the electron 
wave functions corresponding to initial momentum k^ and final momentum kj with respectively outgoing and ingoing 
boundary conditions. The helicity conservation factor Sx^Xi is produced by matrix elements that incorporate the lower 
components. Because the lower components are ±1 times the upper components, in what follows one needs only the 
upper components of the wave function. 

In the plane- wave impulse approximation (PWIA), Coulomb distortion of the electron waves is neglected and the 
electron is described by a plane wave. As a result the integration over r produces (5^^-'(q — Q) with Q = ki — kf. We 
get 

kaPWIA X ^PWIAJNtiiQ, P) 

M ^dxfX, ^2 , (12) 

where h^^j^ denotes known helicity dependent factors (s^]. The PWIA cross section may be expressed in terms of 
longitudinal and transverse response functions, Rl and Rt 

^nJdE, ^ '^Mottj^J i?L + i?T j, (13) 

where 



CTMott = , (14) 



and 

1 , 2Q^ 2^e 



(15) 



With Coulomb corrections included, the longitudinal matrix element of interest must take a gauge invariant form. 
This requires that the electron current must be conserved in the sense that 

d\^i-^\v)(uf, - q . je)e-^'i-*i+)(r) = 0, (16) 
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and that the nuclear current should separately be conserved, 

J d'r^lr^*{r)e'^ '' {iuJ% - q • Jjv) Hr) =0, (17) 

where q is the photon three momentum. With Coulomb distorted waves, the photon momentum q differs from the 
electron's momentum transfer Q = — and the longitudinal current is defined with respect to the direction of the 
photon that is exchanged, not with respect to the difference of asymptotic electron momenta. 
In view of current conservation, the longitudinal current matrix element can be simplified to 

feA - ii ■ m ■ ja^) = (i - ^) ■ (18) 

Using Eq. (O, we get for the longitudinal contribution to A4 



M,= = S,,,^ I d'r I ^e'(Q-<>)re^x(r);^^(-)(r)//^(+)(r)/,0(r)[^i_^ jA'(q,p), (19) 



where Eq. ((TS]) has been used to include the components of je and 3n that are parallel to q. The electron's momentum 
transfer is Q = k^ — k/ and x = X/ ''(i") + Xi^\^) includes the phases of initial and final electron states. Note that 

Xi^\ and 7}^^ = 71^' + iS'l^'^ are obtained from Eq. ([4]) with the z-axis parallel to initial momentum k^. In 

passing, we note that the Glauber approximation is obtained when the eikonal phases for initial and final states are 
evaluated using for each a z-axis parallel to the average momentum, ^(k^ + kj), and only the leading-order phases, 

Xo^^ and Xo~^ ^.re retained. This approximation omits the focusing factors. 

The longitudinal response function is obtained by dividing the cross section integrated over the angles of the 
knocked-out nucleon by the Mott cross section, 

^ / dil„- -r 

(JMottQ^ J (27r)5 



Rl = / dnj,-—^\MLVE}pEp, (20) 



where Ml in the longitudinal amplitude of Eq. (|19p . The full calculation thus involves a six-dimensional integration in 
order to obtain the amplitude Ml- Two more integrations over the angles of the knocked-out nucleon are required in 
order to obtain the response function. Results based on the eight-dimensional integration are called "full calculations" 
in the following sections. The transverse transition matrix element, which will not be considered further in this work, 
is simply the difference of Eqs. ([TO)) and (fTO|) . 

In the actual calculations we mostly use a very simple model for the nuclear current 

where '0(k) is a gaussian wave function for a bound nucleon, 

V'(k) = (2^/32)3/4e-/3^fcV4^ (22) 

normalized such that / d^fc|-0(k)p/(27r)^ — 1. 

This simple gaussian model is used because the Coulomb corrections should depend mainly on the electron wave 
functions. In order to get some idea how a more realistic model of nuclear structure would affect the results we 
considered also shell-model wave functions for ^^Fe and "^^^Pb nuclei. In those cases, the Coulomb potential was 
calculated based on the empirical charge densities of Ref. ^38j] and the range parameter R of the Coulomb potential was 
determined so that the average Coulomb potential matched the empirical one in the sense that J d^ r P(ixpt(T')Vc{r) = 
J d^ r pexpt{f)Vexpt(i') ■ See Table|T]for the parameters used. For the shell-model wave functions, the gaussian parameter 
/3 was selected such that the charge radius of the nuclei agreed with the empirical charge radius and when the higher 
orbitals are included, they are assumed to be described by harmonic-oscillator wave functions in coordinate space as 
follows, 

i^nUv) = NYu^{Slr)r\Fi{-{n^l)/2,l + i/2,rV2/l3)e-^'''^'^\ (23) 
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with normalization constants N determined by J (Prlip^r)]"^ = 1. Furthermore, Yim are the well known spherical 
harmonics and iFi the confluent hypergeometric functions. For the multi-shell studies we in addition use a more 
realistic nuclear current of the form 

= K'^Mp) [^F, + ^^F^cT^^'q,] u{p - q) ^(q - p), (24) 

where Fi{Q^) and F2{Q^) are nucleon form factors, n is the anomalous magnetic moment and K = ]VP / (EpEp^q) is 
a normalization factor arising from the spinors. For the form factors F^ standard dipole ones are taken. 

The nuclear current Eq. (j2ip is based upon the current operator of a scalar nucleon with initial and final momenta 

pl^iEp, p), 

- p-q), (25) 

where uj and q are the photon's energy and momentum. Because of energy conservation, Ep = M + lu — B, where 
B « .OOSGeV is a typical binding energy of a nucleon. For the PWIA response function the angular integration in 
Eq. ([20]) can easily be done. We find for the gaussian model 

yPWIAf.. r\\ 1 {2Ep-UjY (3 f -0^(\Q\-pf/2 -/j^(|Q|+p)V2 



^Piv/A jg normalized so that at fixed Q, J du!RL{Q,uj) ~ 1. 



IV. EFFECTIVE-MOMENTUM APPROXIMATIONS 



Let us consider the electron current matrix element for emission of a photon using spinors corresponding to initial 
and final helicity Xi and A/. Using Eq. ^ and a similar relation for the incoming electron wave function, the electron 
current can be rewritten as 




Using Eq. ([27]) the quasi-elastic transition matrix element ()10p takes the form 

M=S,^,Jd^r J ^e^(Q-'i)V'^W/;^(-)(r)/fW(r)/.^(r)(^i^^ J^^(q,p), (28) 

with 

We may now use a stationary-phase-like argument to calculate Eq. ([^5]) . We see that for large Q the integrand of ([^5]) 
has a rapidly changing phase except when Q — q + Vx(r) = 0. So we expect that the dominant contribution in the 
integrals comes from around this point. Expanding the photon propagator around q = Qe// = Q + Vx(r ~ 0) we 
can explicitly factor the photon propagator out of the integral over q. Hence we expect that Eq. can in a good 
approximation for large Q be determined by 

-M^^^^ = ^A.A. (^^) / ^e^(^-<>)-e^'^W/f (-)(r)/f (+Hr)/.^(r)J,f (r,p) (29) 

with 

Eq. (j29p has obviously the form of an effective momentum approximation except that the full r-dependence of the 
eikonal phase and the focusing factors is retained. It is considerably simpler to calculate than the full six dimensional 
integral for the quasi-elastic matrix element A4. 
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A further approximation can be made by approximating the eikonal phase by x ~ x(0) + r ' ^x(O) but keeping the 
r-dependence of the focusing factors, as foUows, 

M^'^^^- = 5,,,^ (^)^^"^°^ / (02e^(Q-"'>) v;^-^(r)/f (+)(r)/.^f(r)J^^(r,p) (30) 

This is called the EMAr' approximation. Finally one may take both the eikonal phase and focusing factors at the 
central value r = 0. In so doing we get the often-used ema approximation. This approximation usually is based 
on expanding the eikonal phase in a Taylor's series about r = and keeping the first two terms. Moreover, the 
focusing factors are approximated by their values at r = and the helicity matrix elements are approximated by the 
plane- wave values. Integration over r then gives 5^^^(q — Qe//), so the longitudinal amplitude simphfies to the PWIA 
form 

MT" = 2n6x,x,e^^^"^h%^jMQeff,p)^ '-^ —. (31) 

Combining the e^^^ ' « 1 — Vc{0) / {2Ei) factor of the eikonal correction with the (1 — Vc/i?2i)^^^ yields a focusing 
factor Z/'^^'' « 1 — Vc/Ei in the Dirac wave function, thus reproducing at r = the expected factor 1 — Vc(0)/Ei that 
has been derived by Yennie, Boos and Ravenhall |24| based on a WKB analysis of the Dirac-Coulomb wave function. 
A similar result holds for the final-state focusing factor, ff^~\ which is approximately 1 — Vc/Ef. Thus, the overall 
focusing effect in the matrix element is approximately equal to (1 — Vc(0)/£'/)(l — V{0)/Ei). 

The effective momentum involves the gradient of the eikonal phase shift x — x') ' + xl^' the origin. Because of 
cylindrical symmetry of x^^'' about the direction ki, Vx^^'' at the origin is nonzero only along the direction ki, and 
similarly Vxy at the origin is nonzero only along the direction kf. With Vi = Vf ~ 1, we find the same result as 
Traini, 

Qe// = ki ki — 6k — kf kf — Sk , (32) 

where Sk = Vc(0). It is correct up to first order in the eikonal expansion because the contribution from the gradient 
of eikonal correction xi vanishes at the origin. 

As shown by Rosenfelder [2^ and Traini [2^ . there are significant cancellations in the Coulomb corrections when 
response functions are evaluated in this effective-momentum approximation (ema) using the approximate focusing 

factors, f[""^^ « 1 — Vc{0)/Ei and "* ~ 1 — Vc{0)/Ef. Coulomb effects in the focusing factors and the effective 
photon propagator cancel if one considers the photon propagator of the transverse amplitude, which is l/[Qe ' '^^ ^ 



l/[4[fc,-y,(0)][fc/-l^e(0)]szn2i 



1-6., 



Qlff--^' ~ Q'- ^^^^ 

These factors produce the same result as in the plane- wave case, Eq. (|12p . but the momentum transfer argument in 
the nuclear structure function is shifted. 

Evidence has been presented that the momentum shift as predicted by the ema is too large and that a smaller value 
should be taken for the Coulomb potential at the origin. This is done in view of the plausible classical argument that 
the Coulomb potential which is felt by the electron is not the central value of the potential, but rather is the average 
potential along the electron trajectory. Based on this argument the momentum shift 5k in the eikonal wave function 
is weakened by a factor fema 

Sk = femaVM, (34) 

where fema is determined by fitting the experimental quasi-elastic peak value. In practice one finds a reduction factor 
of typically fema « 0.7 to 0.8. 

It should be noted that although the Qeff is modified by the factor fema, in the actual analysis one assumes 
that the cancellation ([55)1 still holds. In general, this cancellation is clearly expected not to be complete. Using a 
gauge-invariant response function as obtained from Eq. (jl9[) leads to deviations, which are of the order of one percent. 
Another source of deviation is a more precise treatment of the r-dependence of the focussing factor as is done in the 
EMAr approximation, given by Eq. (j29p . The break down of the cancellation in the effective-momentum approach is 
reflected in allowing for an additional overall normalization in the response function as given in Eq. psp . 
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TABLE I; Parameters used in calculations: /3 is the harmonic oscillator parameter; Vb and R are the Coulomb potential 
parameters. The gaussian model refers to Eq. (|22p . 



Nucleus 


/3(fm) 


Vo (GeV) 


R (fm) 


gaussian 


2.0 


0.0273 


2.0 


208 


3.564 


0.0256 


7.10 




2.854 


0.0124 


3.97 



V. RESULTS 

In this paper, we have described the eikonal expansion for relativistic wave functions in the presence of a Coulomb 
potential based on the Dirac equation. In the considered expansion, focusing factors are obtained in a systematic 
manner by use of the eikonal expansion. Although focusing factors take somewhat different forms u sing the Klein- 
Gordon wave function, equivalent results are found for the current matrix elements for the two cases [30| • 

Calculations of the longitudinal response function are performed for four cases: PWIA, ema , EMAr and the full 
calculation using distorted waves based on the Dirac equation. Eikonal phases are evaluated through second order, 
i.e., X = Xo + Xi + X2 and lu — uji + UJ2- It should be noted however, that the expansion converges rapidly for the 
parameters and energies used and results based on xo + Xi and cui differ by about 0.3% at the quasi-elastic peak. 
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FIG. 2: Longitudinal response function versus the electron's energy loss, u, calculated using the gaussian model for e~ scattering 
at i5 = 500MeV and 0^ = 60°. Dotted line shows PWIA, solid line shows ema based on fema ~ 0.7 and x's show full calculations 
based on Eq. ((19}. 

In the figures we have used the gaussian model of Tableland have considered electron and positron scattering. 
Calculations have omitted final-state interactions of the knocked-out nucleon. Figure[2] shows the longitudinal response 
function for 500 MeV electrons and electron scattering angle 9^ = 60°. Here the full calculations are plotted as x's, 
while the ema calculation based on 6k = 0.7 14(0) is shown by the solid line. The PWIA prediction is shown by the 
dotted line. Figure [3] shows similar results for the longitudinal response function for e"*" scattering at 540 MeV using 

fema — 0.7. 

In general the ema is seen in Figures [5] and [3] to produce a significant shift oi Rl away from the PWIA result 
and towards the full calculation of R^. In both cases, the full calculations are reproduced quite well by the ema 
using fema — 0.7. This reconfirms the findings of Refs. [13, HI] that a smaller value of the 5k than Vc(0) produces 
better agreement with the full results. From the effective-momentum approximation results we see that the shift in 
momentum due to the Coulomb distortion is predicted to be opposite in e+ to that of e— scattering. There is indeed 
good agreement between the response functions for e— and e+ scattering at the energies that make Qe// close to the 
same for both. 

As discussed in Ref. [30| , it is possible to fit the response functions more precisely if the momentum shift is allowed 
to be a function of energy loss, lo. For the gaussian model, it is found that near the quasi-free peak the momentum 
shift is well described by fema — 0.7. The momentum shift can be significantly larger in magnitude, corresponding 
to fema > 1, whcn LO is significantly away from the value at the quasi-free peak. The reason is that the Coulomb 
distortions tend to alter the shape of the response functions away from the quasi-free peak. However, simply using a 
constant /ema does not incur large errors. The response integrated over oj, as in the Coulomb sum rule, is expected 
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FIG. 3: Longitudinal response lunction versus the positron's energy loss, u, calculated using the gaussian model for 
scattering at _E = ^AQMeV and 6e = 60°. The dotted line shows PWIA, solid line shows ema using fema = 0.7 and x's show 
full calculations based on Eq. (fT9)l. 



TABLE IL The full calculation of the response function for the gaussian model at various uj values for e scattering at 500 
MeV, 9e ~ 60°, together with the effective-momentum approximation results. 





Full 


ema 


EMAr 


EMAr' 


0.060 


1.616 


1.678 


1.692 


1.756 


0.080 


4.548 


4.555 


4.603 


4.655 


0.100 


7.626 


7.644 


7.672 


7.717 


0.120 


9.122 


9.083 


9.066 


9.122 


0.140 


8.236 


8.292 


8.280 


8.326 


0.160 


6.124 


6.144 


6.191 


6.194 


0.180 


3.854 


3.847 


3.955 


3.909 


0.200 


2.124 


2.100 


2.231 


2.160 


0.220 


1.055 


1.025 


1.143 


1.075 



to be accurate within one or two percent. 

In Table |TT] we show the numerical results for scattering ai E = 500MeT^ and 9^. = 60" of the full calculation 
together with the various effective-momentum approximations, using the gaussian model. From this we see that 
there is a close agreement between EMAr and the approximation EMAr', obtained from EMAr on replacing the 
eikonal phase by x{f) ~ fema r • Vx(0). This illustrates that the r-dependence of the phase shift can indeed well be 
approximated by the linearized form. Moreover, both ema and EMAr are in good agreement with the full result. It 
should be noted that the assumed complete cancellation of the focus factor (j33[) in the ema approximation, which is 
found to hold in this case, may be accidental. 

The results for the simple gaussian model suggest that the ema can reproduce the results of the more elaborate 
EMAr analysis quite well. In order to test this proposition for a more realistic model of the quasi-free scattering, 
calculations have also been made for the ^^Fe and ^'^^Pb nuclei using shell-model wave functions and a Coulomb 
potential that is based on the empirically determined charge density. The Dirac nucleon current Eq. ([M)) is used. 
Figures [3] and [5] show the results. 

In these figures, we show the EMAr calculations as solid lines and two-parameter fits to them using the ema-fitting 



TABLE III: Parameters used to fit the EMAr response functions using Eq. (|35p. 



Nucleus 


E 


<? 


5k 


fema 


A 


58 
208 


0.5 
0.5 


0.55 
0.55 


-8.8 
-21.0 


0.71 
0.82 


0.99 
0.98 
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0) (GeV) 

FIG. 4: EMAr longitudinal response function for ^'^Fe at Ei — 500 Mev and q = 550 MeV/c (solid line). The corresponding 
PWIA response functions without Coulomb effects included is shown by the dotted line. A fit of the EMAr response function 
using Eq. (|35|l is shown by the x symbols and values of the fitting parameters are given in Table Hill 




CO (GeV) 

FIG. 5: EMAr longitudinal response function for ^"^Pb at Ei = 500 Mev and q = 550 MeV/c (solid line). The corresponding 
PWIA response function without Coulomb effects included is shown by the dotted line. A fit of the EMAr response function 
using Eq. (|35|l is shown by the x symbols and the parameters of the fit are given in Table Hill 



formula of Eq. ([35|) as x 's. 

Rl{Q,uj) - A^"^^ i?r^^^(Qe//,c.) (35) 

The obtained parameters are given in Table IIIII The main fitting parameter is the value of fema that is used to 
determine the momentum shift in Qe//. Using the shell-model wave functions, we find that response functions are fit 
by fema = 0.71 for ^^Fe and by fema = 0.82 for ^^^Pb. Because of small distortions of the shape of the EMAr response 
function relative to the shape of the PWIA response function, a minor change of normalization is used also, as given 
by the A parameter. The Coulomb effects are larger for the ^°^P6 nucleus because of the larger Coulomb potential, 
however the results can be fit using the ema formula with an appropriate value of fema ■ The fact that similar values 
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0.1 
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0.3 
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CO GeV 



FIG. 6: Solid lines show the partial EMAr longitudinal response function for the Ip, Id and Ig shells of '^"^Pb at = 500 
Mev and q = 550 MeV/c. The corresponding ema response functions based on fema = 0.82 are shown as dash lines and the 
PWIA response functions are shown as dotted lines. The upper group is the Ig shell, the middle group is the Id shell and the 
lower group is the Ip shell. 

of fema are found for the gaussian model, for ^^Fe and for "^^^Ph demonstrates that the Coulomb corrections are not 
very sensitive to the nuclear model. 

Figure [6] shows partial response functions for individual shells of ^^^Pb at 500 MeV electron energy. Three shells, 
the Ip, Id and Ig shells of ^^^Ph, are shown based on three calculations: the EMAr, ema using fema = 0.82 and 
PWIA. The results for individual shells show that the ema results are close to the EMAr response functions that take 
into account the r-dependence of eikonal and shell wave functions. Differences are somewhat larger than for response 
functions summed over all shells. This is expected because the radial wave functions for the shells with € > are 
suppressed near the origin where the Coulomb potential is largest. The ema uses an average value for the Coulomb 
potential that works well for the sum over all shells and is less accurate for individual shells. 



Within the eikonal approach we have studied how well the Coulomb distortion effects are described by effective- 
momentum approximations in electron scattering on nuclei at intermediate energies. This is expected to be a reliable 
description at increasing high energy. We have shown using a systematic eikonal expansion that already in the 
few hundred MeV region the convergence of the expansion is indeed very fast and that the leading orders of the 
eikonal phase and focus factors in the electron wave function are sufficient to describe the Coulomb distortion in an 
accurate way. Moreover, the effective-momentum approximations are found to agree well with the full eikonal based 
calculations. 

In this paper we have focused on the longitudinal response function. The transverse response contribution has also 
been calculated together with the spin-dependent terms occurring in the eikonal wave function [sj ■ From the present 
study we find strong support for the conjecture that the effective-momentum approach can be used as the basis for 
analysis of the inclusive experimental data. In particular, the ema approximation as used in the actual analysis of 
the experimental data is found in our model studies to do well. 

A more precise form of the effective-momentum approximation would be useful for removing Coulomb corrections 
from experimental data in a straightforward manner and has been suggested in Ref. (30| . In order to have a precise 
result, one could determine appropriate values of the momentum-shift function Sk(ki, uj, 9e) from which the appropriate 
Qe// may be calculated as in Eq. ([5^ . 



VI. CONCLUDING REMARKS 
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In the ema approximation a weakening factor fema for the momentum shift is usually introduced, which takes into 
account in a phenomenological way the trajectory dependence of the phase shift. We have studied in this paper the 
EMAr approximation, which includes explicitly this r-dependence in the eikonal shift and focus factor without the 
introduction of fema- Up to a possibly small overall normalization constant correction, this is found to be in good 
agreement with the full calculation 

Our results for response functions omit final state interactions of the knocked-out nuclcon. Their inclusion would 
affect the shape of the response functions but not the value of the Coulomb sum rule, which involves a sum over a 
complete set of states. The Coulomb corrections found in this work are not large enough to explain the differences 
that have been reported for the Coulomb sum rule by different experimental groups. jol. [ill [l3| 
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